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Abstract. The inverse spectral problem is investigated for the matrix Sturm-Liouville 
equation on a finite interval. Properties of spectral characteristics are provided, a con- 
structive procedure for the solution of the inverse problem along with necessary and 
sufficient conditions for its solvability is obtained. 
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1.1. In this paper, the inverse spectral problem is investigated for the matrix Sturm- 
Liouville equation. Inverse spectral problems are to recover operators from their spectral 
characteristics. 

The scalar case has been studied fairly completely (see [1-3]). The matrix case is more 
difficult for investigating. Different statements of inverse spectral problems for the matrix 
case were given in [4] , [5] and [6] with corresponding uniqueness theorems. A constructive 
solution procedure was provided in [7] , but for the special case of the simple spectrum only. 
Necessary and sufficient conditions were obtained in [8] for the case when the spectrum 
^ ■ is asymptotically simple, that is an important restriction. Moreover, the method used by 
On , the authors of [8] does not give a reconstruction procedure. We also note that necessary 
and sufficient conditions on spectral data were given in [9] for Sturm-Liouville operators 
'vh I with matrix- valued potentials in the Sobolev space . This class of potentials differs 

00 ■ from one considered in this paper. 

In this paper, we study the self-adjoint matrix Sturm-Liouville operator in the general 
case, without any special restrictions on the spectrum. Properties of spectral charac- 
teristics are investigated, and necessary and sufficient conditions are obtained for the 
solvability of the inverse problem. We provide a constructive procedure for the solution of 
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' the inverse problem in the general case, that is a generalization of the algorithm from [7]. 
^ For solving the inverse problem we develop the ideas of the method of spectral mappings 

[3]. 

1.2. Consider the boundary value problem L = L{Q{x), h, H) for the matrix Sturm- 
Liouville equation: 

iY := -Y" + Q{x)Y = XY, a;e(0,7r), (1) 

U{Y) := Y\0) - hY{0) = 0, V{Y) := Y\n) + HY{tt) = 0. (2) 

Here Y{x) = [yk{x)]i.^Y^ is a column vector, A is the spectral parameter, and 
Q{x) = [Qjk{x)]j f.^i^ , where Qjk{x) G L2(0,7r) are complex-valued functions. We will 
subsequently refer to the matrix Q{x) as the potential. The boundary conditions are 
given by the matrices h = [hjk]j^k=l~m > ^ = [Hjk]j^k=l~m > where hjk and Hjk are 
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complex numbers. In this paper we study the self-adjoint case, when Q — Q* , h — h* , 
H = H* . 

Let if{x,X) and S{x, \) be matrix-solutions of equation (1) under the initial condi- 
tions 

cpiO, A) = Im, ^'{0, A) = h, S{0, A) = 0^, S'{0, A) = 7^. 

where 1^ is the identity m x m matrix, 0^ is the zero m x m matrix. 

The function A(A) := det[K(v9)] is called the characteristic function of the boundary 
value problem L . The zeros of the entire function A (A) coincide with the eigenvalues of 
L (counting with their multiplicities, see Lemma 3), and they are real in the self-adjoint 
case. 

Let u — u!* be some mxm matrix. We will write L{Q{x),h,H) & A{ijj) , if the prob- 
lem L has a potential from ^2(0,7?) and h+H+^ Q{x) dx = uj . Without loss of gen- 
erality we may assume that L e A{uj) , where cu & D — {cu: cu — diag{a;i, . . . , oum}, < 

. . . < UJm} ■ 

One can achieve this condition applying the standard unitary transform. 
In order to formulate the main result we need the following lemmas that will be proved 
in Section 2. 

Lemma 1. Let L e A{uj) , ou & D . The boundary value problem L has a countable 
set of eigenvalues {Xng}n>o,q=T^ > ^'^^ 

Pnq= ^/\^=n+ — + —, {Knq}n>0 ^ I2, q = '^,m. (3) 

irn n 



Let A) = [^jk{x, A)]j j.^]^ be a matrix-solution of equation (1) under the bound- 
ary conditions f/($) = Im , = 0™. • We call A) the Weyl solution for L . Put 
M(A) := $(0, A) . The matrix M(A) = [Mjki\)]j^k=l-^ is called the Weyl matrix for L . 
The notion of the Weyl matrix is a generalization of the notion of the Weyl function ( m - 
function) for the scalar case (see [1], [3]). The Weyl functions and their generahzations 
often appear in applications and in pure mathematical problems, and they are natural 
spectral characteristics in the inverse problem theory for various classes of differential 
operators. 

Using the definition for M(A) one can easily check that 

M(A) = -{V{^)r'V{S). (4) 

It follows from representation (4) that the matrix-function M(A) is meromorphic in A 
with simple poles in the eigenvalues {Xnq} of L (see Lemma 4). 
Denote 

anq := Res M(A). 

X=Xnq 

The data A := {Xnq, <^nq}n>o,q=l~m called the spectral data of the problem L . 
Let {Xn^qi^}k>o be all the distinct eigenvalues from the collection {Xnq}n>o,q=l~m ■ 

^'nkqk '■= ^^klk^ ^ ^ 0' ^'nq = ^rn, {u, q) ^ {{Uk, qk)}k>0- 

Denote 

1 = mi < . . . < rup+i —m + 1, 
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where p is the number of different values among {oUq}g^i^ . Let 



rus+i-l 
q=ma 



Lemma 2. Let L e A{uj) , uj E D . Then the following relation holds 

+ {||4^)||}n>oe/2,5 = T;^, (5) 



TT n 



where 



r(s) ^ rAsh r{s) ^ J 1> rns<j = k< m^+i - 1, 

[ jfc Jj,fc=i,m> jfe I 0, otherwise, 



and ||.|| a matrix norm: \\a\\ = maxj^ajk 
Consider the following inverse problem. 

Inverse Problem 1. Given the spectral data A , construct Q , h and H . 

We will write {A„g, oinq}n>o,q=T;^ e Sp , if for X^g = A^; we always have a^g = aki . 
The main result of this paper is 

Theorem 1. Let uj E D . For data {Xnq, (^nq}n>o,q=l~m ^ spectral data 

for a certain problem L e ^(<^) it is necessary and sufficient to satisfy the following 
conditions. 

1) The asymptotics (3) and (5) are valid. 

2) All Xng are real, ang — {(Xnq)* , ocng > for all n >0 , q — l,m , and the ranks 
of the matrices ang coincide with the multiplicities of Xnq ■ 

3) For any row vector 7(A) that is entire in X , and that satisfy the estimate 

7(A) = 0(exp(|/m\/A|7r)), |A| ^ 00, 



if l{.^nq)otnq = for all u > Q , Q = l,m , then 7(A) = . 

We prove necessity of the conditions of Theorem 1 in Section 2 and sufficiency in 
Section 4. In Section 3 the constructive procedure is provided for the solution of Inverse 
Problem 1. 

2. Necessity. 

2.1. Let us study some properties of the spectral data. 

Lemma 3. The zeroes of the characteristic function A(A) coincide with the eigen- 
values of the boundary value problem L . The multiplicity of each zero Aq of the function 
A(A) equals to the multiplicity of the corresponding eigenvalue (by the multiplicity of the 
eigenvalue we mean the number of the corresponding linearly independent vector eigen- 
functions). 

Proof. 1) Let Aq be an eigenvalue of L , and let be an eigenfunction correspond- 
ing to Ao . Let us show that Y^{x) = cp{x, Xo)Y^{0) . Clearly, Y^{0) = (^(0, A)y°(0) . 
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It follows from U{Y^) = that ^^'(O) = hY^{0) = ip{0,X)Y^(0) . Thus, Y^x) and 
ip{x, Xo)Y^{0) arc the solutions for the same initial problem for the equation (1). Conse- 
quently, they are equal. 

2) Let us have exactly k hnearly independent eigenfunctions Y^ , , . . . , Y'' 
corresponding to the eigenvalue Aq . Choose the invertible m x m matrix C such that 
the first k columns of ip{x,\o)C coincide with the eigenfunctions. Consider Y{x,\): = 
(p{x,X)C , Y{x,\) = \Yq{x,X)\q^j^^, Yq^x.Xo) = Yi{x),q = T/k . Clearly, that the 
zeros of Ai(A) :— det V{Y) — det V{(p) ■ det C coincide with the zeros of A(A) counting 
with their multiplicities. Note that A = Aq is a zero of each of the columns V{Yi) , 

V{Yk) . Hence, if Aq is the zero of the determinants Ai(A) and A(A) with the 
multiplicity p , than p > k . 

3) Suppose that p > k . Rewrite Ai(A) in the form 



As (A) = det 



Ai(A) = (A-Ao)'=A2(A), 
V{Y,) V{Yk) 



,y(n+i),...,y(y^) 



A — Ao A — Ao 

In view our supposition, we have A2(Ao) = , i. e. there exist not all zero coefficients 
aq , q — l,m zero exist such that 

^^^ dy(n(..Ao)) ^ ^ a,V(Y,(x,X«)) = 0. (6) 

q=l q=k+l 



If Q!g = for q = l,k , then the function 

m 

Y+{x,X):= J2 

q=k+l 

for A = Aq is an eigenfunction corresponding to Aq that is linearly independent with 
y , q = l,k . Since the eigenvalue Aq has exactly k corresponding eigenfunctions, we 
arrive at a contradiction. 

Otherwise we consider the function 

k m 

Y+{x, A) := Yl "^M^^ A) + (A - Aq) ^ aqYq{x, A). 

q=l q=k+l 



It is easy to check that 

£{Y+(x, A)) = XY+{x, A), i [^Y+(x, X) j = A^y+(x, A) + Y+(x, A), 

C/(y+) = U (^^y+) = 0, V{Y+{x, Ao)) = 0. 
Relation (6) is equivalent to the following one 
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Thus, we obtain that Y^(x,Xo) is an eigenfunction, and -^Y^(x, Xq) is a so-called 
associated function (see [10]) corresponding to Aq . If we show that the considered Sturm- 
Liouville operator does not have associated functions, we will also arrive at a contradiction 
with A2(Ao) 7^ , and finally, prove that k — p . 

4) Let us prove that the self-adjoint operator given by (1), (2) does not have associated 
functions. Let Aq be an eigenvalue of L , and let and Y^ be a corresponding 
eigenfunction and an associated function respectively, i. e. both F° and Y^ satisfy (2) 
and 

(^-Ao)y° = o, (£-Ao)y' = y°. 

This yields 

{{£ - XorY\Y') ^ 0, 

for the scalar product defined by 



PIT 

{Y,Z) := / Y*{x)Z{x)dx. 
Jo 



In case of the sclf-adjoint operator, we have {iY,Z) — [Y,£Z) for any Y and Z satis- 
fying (2), and the eigenvalue Aq is real. Therefore, 

((£-Ao)yS(^-Ao)y^) = (nn = o, 

and = . Recall that Y^ is the eigenfunction, and get a contradiction. □ 

Lemma 4. All poles of the Weyl matrix M(A) are simple, and the ranks of the 
residue-matrices coincide with the multiplicities of the corresponding eigenvalues of L . 

Proof. Let Aq be an eigenvalue of L with a multiplicity k , and let Yi , Y2 , . . . , Y^ 
be linearly independent vector eigenfunctions corresponding to Aq . Following the proof 
of Lemma 3, we introduce the invertible matrix C — [Ci, . . . ,Cm\ such that Yq{x) ~ 
(f{x,Xo)Cq , q = l. k . Consider the vector-function Y(x,X) — (p{x,X)C . Clearly, that 
{V{(p))-^ = C{V{Y))-^ . Write V{Y{x,X)) in the form 

V{Y{x, A)) = [(A - Ao)l^i(A), . . . , (A - Xo)Wk{X), Wk+i{X), W^(A)], 

where 

V{Y,{x,X)) 



A — Aq 



Wq{X) = V{Yq{X)), q^k + l,m. 
Clearly, Wq{X) are entire functions, and 

detW{X) = det[W^i(A), . . . , Wm{X)] ^ 

for A from a sufficiently small neighborhood of Aq (otherwise the multiplicity of the 
eigenvalue Aq is greater than k \ It is easy to show that 

k 



det V{y{x, A)) = (A - Ao)^ det W{X), 



{y{Y{x,X)))-' 



^i(A) Xfe(A) 



A — Aq a — Aq 



(A),...,X„(A) 
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where Xq{\) are analytic in a sufficiently small neighborhood of Aq (the superscript t 
stands for transposition). Using (4) we get 

ao = Res M(A) = - Res{V(ip(x, X)))-^V(S(x, A)) 

A=Ao A=Ao 



= - Res C 

A=Ao 



^ ^ ^ ^ -^fc+i(A), . . . , Xto(A) 



V{Six,X)) 



A — Ao A — Ao 

= -C [Xi(Ao), . . . , Xk{Xo), 0, . . . , 0]* ViSix, Ao)) = -XV{S{x, Ao)). 

Therefore, the poles of the Weyl matrix are simple, and rank ao < k . 
Let us prove the reverse inequality. Note that 

Res{V{^{x, X)))-'V{^{x, A)) = 0^ = XV{^{x, Aq)). 

A=Ao 

Let ip{x,Xo) be a solution of equation (1) for A = Ao under the initial condition V{ip) — 

X* . Since columns of the matrices (f{x, Ao) and S{x, Ao) form a fundamental system 
of solutions of equation (1), we have 

iIj{x, Aq) = ip{x, Xo)A + S{x, Xo)B, 

XX* = XViiPix, Ao)) = XV{(p{x, Xo))A + XV{S{x, Xo))B = -aoB. 

On the one hand, since det W ^ , the vectors Xg(Ao) are linearly independent, there- 
fore, rank XX* = k . On the other hand, rank agB < rank ag . Thus, we conclude that 
rank qq> k . □ 

Lemma 5. Let Xq , Ai he eigenvalues of L , Xq ^ Xi , and ctj = Res M(A) , 

i = 0, 1 . The following relations hold 



-K 

j ^*{x,Xq)v{x,Xo) dxao^al, 







TT 

J ip*{x,Xo)ip{x,Xi)dxai^Or, 





In particular, according to the first relation, 

OiQ — a.Q> 0. 

Proof. Denote 

tZ := -Z" + ZQ{x), V*{Z) := Z'(7r) + Z{t:)H, (Z, Y) := Z'Y - ZY' , 
where Z = [Zk]\^YTa ^ vector ( t is the sign for the transposition). Then 

= V*{Z)Y{7r) - Z{7r)V{Y). 
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If Y(x,X) and Z{x,fj,) satisfy the equations iY(x, \) = XY(x, X) and i*Z(x,fj,) = 
IJ,Z{x, ii) , respectively, then -^{Z, Y) = {X — ii)ZY . In particular, if X = /i , then {Z, Y) 
does not depend on x . 

Since Aq is real, (p*{x,Xq) satisfies the equation i*Z — XqZ . Hence, 



/ (p*{x, Xo)(p{x, Xq) dx — lim 
Jo ^^^0 



lim 



{^*{x,Xo),^{x, A))|g 

A— >-Ao A — Aq 

V*{ip*{x, Xo)Mx, X) - v*{x, Xo)V{ip{x, A)) 



A— s-Ao A — Aq 

It follows from (4) and Lemma 4 that 

V{ip{x, Xo))ao = - lim (A - Xo)V{ip{x, X)){V{if{x, X)))-'V{S{x, A)) = 

A->-Ao 

Analogously aQV*{ip*{x, Xq)) = Om ■ Consequently, we calculate 

*/*/ \ \ / \ \7 **/ \ ''{^{■^■>X)) 

Q!o / (fi {x,Xo)(p{x,XQ)dxaQ = aQifi (tt, Aq) hm — 

J A-^Ao A — Aq 



X lim (A - Xo){V{if{x, X))y'V{S{x, A)) = a*(^*(7r, Xo)V{S{x, Aq)) 

A^Ao 

= -tto(v'*(^' Ao), S'(a;, Xo))x=Tr = -aQ{Lp*{x, Xo),S{x, Xo))x=o = "o- 

Similarly one can derive the second relation of the lemma. □ 

2.2. In this subsection we obtain asymptotics for the spectral data. 
Denote p := y/X , Re p > , r := Imp , Gg = {p: \p — k\ > S,k = 0,1,2, ...} , 
5 > . By the standard way (see [3, Sec. 1.1]) one can obtain the estimate 

A(A) = (-psinp7r)'"+0(|pr-^exp(m|T|7r)) = (-psinp7r)"+(-psinp7r)"-^0(exp(|r|7r)) 

+ . . . + (— psinp7r)0(exp((m — l)|T|7r)) + 0(exp(m|T|7r)), |p| — >■ oo. (7) 

Proof of Lemma 1. 1) Consider the contour = {A: |A| = + . By virtue 

of (7) 

A(A) = /(A) + ^(A), /(A) = (-psinpTT)™, |^(A)| < Clpp-^ exp(m|r|7r). 

If A e Fat for sufficiently large , we have |/(A)| > |5'(A)| . Then by Rouche's theorem 
the number of zeros of A(A) inside F^v coincide with the number of zeros of /(A) 
(counting with their multiplicities), i.e. it equals (A^ + l)m . Thus, in the circle |A| < 
(A" + 1/2)^ there are exactly (A^ + l)m eigenvalues of L : {Xnq}n=o;N,q=T^ ■ 

Applying Rouche's theorem to the circle 7n((^) = {p- \p — < 6} , we conclude that 
for sufficiently large n there arc exactly m zeros of A(p^) lying inside 7n((^) , namely 
{Pnq}q=T^ . Siuce (5 > is arbitrary, it follows that 

Pnq^ri + Snq, Snq = o(l), U ^ OO. (8) 

Using (7) for p = Pnq , we get 

{-Pnq sin Pnq'^)"^+{-Pnq Sin Pnqn)"^~^0{l) + . . .+{-pnq Sm p„qn)0{l)+0{l) = 0, 71 ^ OO. 
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Denote Snq ■— |pngSinp„q7r| , and rewrite the obtained estimate in the form 

< Co + C,s^, + ... + C^-iCr'- (9) 
It follows from (9) that s^q < maxjl, Ylh^o ^k} • Otherwise we arrive at a contradiction: 



rn—l m—1 

n, e 



k=0 k=0 

Hence, \pnq sin pnq7r\ < C . Using (8) we get 

sinp„g7r = sin£„q7rcosn7r = 0(n~^), £nq = 0{rr^), n — >■ oo. 
Together with (8) this yields 

Pnq^ rb^-0{n~^), n ^ oo. 
2) Let us derive the more precise asymptotic formula. One can easily show that 
V{(fi) — — psin pTT • /to + a; cos fnr + k{p), 

where 

k{p) = Qit) cos pi-K -2t) dt + Qexp(|T|7r)^ . 



Consider the linear mappings Zn{p) that map the circles {p: \p — n\ < C/n} (note 
that pnq lie in these circles for a fixed sufficiently large C ) to the circle {z: \z\ < R} : 

p = n-\ —. 

nn 

For \z\ < R we have 

V{ip) = {-iriuj - Zn{p)Im + l^n{Zn{p)))- (10) 

Using the representation for k{p) wc get k„(z) = o(l) , n — >■ oo , uniformly with respect 
to z in the circle {z: \z\ < R} . Moreover, for each sequence {-2°}n>o C {z: \z\ < 
R} we have {||'^n(^n)ll}n>o £ h and ^ ||'*n(^n)ll^ < C, where C is some constant. 

n>0 

Consequently, 

Aip') = ±fiZnip))+gn{Znip)), 

where f{z) = det(a; — zl^) , gn{z) = o(l) , n — )■ oo (uniformly with respect to z £ 
{z: \z\ < R} ), and the choice of sign ± depends only on n . Fix < S < 1/2 min \u}q— 

q,l:uiq^UJi 

iUil and introduce the contours 7g = {-z: |-2 — Wq| = S} . Clearly, the inequality \f{z)\ > 
\gn{z)\ holds on 7q for all sufficiently large n , and by Rouche's theorem two analytic 
functions A{p'^{z)) and f{z) have an equal number of zeros inside 7^ (here Pn is the 
inverse mapping to Zn)- Thus, we have 

^q ^nq 



Pnq^ri-\ \ , = 0(1), n ^ 00, 1,171. 

Tm n 



9 



Substituting this formula into (10) we get 

Since {||'tn(-2ra(prag))||} ^ h , oue can easily prove that {/t„g} Eh- □ 

Proof of Lemma 2. 1) Let M(A) be the Wcyl matrix for the problem L{Q, h, H) , 
such that Q{x) = , h = H = 0,^ . Then = f /(^) , s = T7p . 

Consider the contours = {A: |A — (n^ + = R} , R = - min \(jOq — 0Ji\ . 

Using the residue theorem and taking Lemma 1 into account, we deduce 

^ „ ms+i-l m^+i-l 

— J (M(A) - M(A)) dA = J] - J] = - -/(^), n>n\s = i;^. 



q=ms q=ms 



7A 

One can easily show that Mjk{X) = — , where 

Using this representation, we arrive at 

M,.(A)-Af,.(A)^^^"^^---(^^-^---(^)^(^\ j,k = Y^. (11) 
jfcv ; jkK ) A(A)A(A) ' ' v ^ 

Let us use the mappings Zn introduced in the proof of Lemma 1: 

, ^n(p) 

p^n-\ —. 

nn 

If A G , then < 6i < \zn{p) — uj,j\ for all q = l,m , and |-2„(p) — oo^^^l < 82 ■ 
Hence, the estimate for A (A) obtained in the proof of Lemma 1 is valid: A(A) — 
±f{zn{p)) + 0(1) , A e 'fn^ , n — )■ 00 (uniformly with respect to A ). 
Similarly, we estimate 

A,.(A)^± ^I'f^^ +0(1) forj^k, 

A,fc(A) = o(l) fcyrj^k, 

A e n — >■ 00, j,k — 1, m. 

Convergence of the remainders is uniform with respect to A , the choice of sign ± depends 
only on n . Analogous relations hold for A(A) and Ajfc(A) . 

Substituting these estimates into (11) and taking into account that Ci < \f{zn{p))\ < 
C2 for A e 'fn^ , we arrive at 

Mjk{X)-Mjk{X)^o{l), j,k^T;^, Xejl:\ 
^ J (M(A) - M(A)) dX = 0(1), ai^) = h^^^ + rj^\ r)^^^ = o(l), n ^ oc. 
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2) Below one and the same symbol denotes various matrix sequences such that 

e I2 . Using the standard asymptotics 

X 

^ , ,sinpx f sin p(x — 2t) ^ , . , ^/explrlxX 
X, A = cospx ■ Im + Qi{x) ^ + / '-Q{t) dt + 0[ ' , 

P J \ P J 



IpI — >■ 00, X e [0, tt], 

a; 

where = h + J Q{t) dt , one can easily show that 



TV 

/TT h\i t\j 

^P*{X, \nq)^{x, Xnl) dx = -/^ + — , X^q - Kl = — ■ 
^ Tt Th 



Applying Lemma 5, we get 

(^^m + ^) "ng = "ng, « > 0, q= 1,171. 



Clearly, ||q;„5|| <C, n>0, gf = l,m. Consequently, fa^^ = anq + ^ ■ Similarly we 
derive anqOCni = ^ , rUg <q,l < m^+i -1, ^t^Z, s = l,p. Thus, 



g=ms / q=ms q=ms 

Substitute the result of point 1 into this equality: 

2 \7r / TT n 

Consequently, 77!*^ = if • ^ 

^.5. Proof of Theorem 1 (necessity). The first two conditions are fulfilled by 
Lemmas 1, 2, 4, 5. 

Let 7(A) be a function described in condition 3. Recall that 

V{(p{x,Xnq))anq = 0^. 

Since 

rank V{(p{x, A„q)) + rank anq = m 

and 7(A„g)a„g = , we get 7(A„g) = CnqV{(p{x, X^q)) , i.e. the row 7(A„q) is a linear 
combination of the rows of the matrix V{(p{x, Xnq) (here Cng is a row of coefficients). 
Consider 

f{X)^^{X){V{^{x,X))r\ 
The matrix-function {V{<^{x, X)))~^ has simple poles in A = A„q , therefore, we calculate 

Res /(A) = 7(A„,) Res {V{^{x,X)))-' 
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= C„, lim V{^{x,X)) lim (A - A„,)(V^((^(x, A)))"^ - 0. 

Hence, /(A) is entire. It is easy to show that 

||(V^((^(a;,A)))-i <C5|p|-^exp(-|r|7r), p e Gg, 

where Gs — {p:\p — k\ > 5,k — 0, 1, 2, . . .} , 5 > . Prom this we conclude that 
||/(A)|| < in . By the maximum principle this estimate is valid in the whole 
A -plane. Using Liouville's theorem, we obtain /(A) = . Consequently, 7(A) = . □ 

Note that in the scalar case condition 3 follows from the first two conditions of Theo- 
rem 1. Indeed, in the scalar case, we have 7(A„)q;„ = , n > , where «„ are positive 
real numbers. Hence, 7(A„) = . Having the spectrum {A„}„>o we can construct the 
characteristic function (see [3, Theorem 1.1.4]): 



A(A) = 7r(A - Ao) n 

n=l 



11? 



and using asymptotics (3) for the eigenvalues we get the estimate 

||A(A)||>C5|p|exp(|T|7r), p e G^. 

Then we introduce /(A) = and follow the proof of necessity in Theorem 1. 

In the general case, condition 3 is essential and cannot be omitted, that is shown by 
the following example. 



Example 1. Let m — 2, Aqi 7^ A02 , A„i = A„2 — 'n^, ^t- > 1 , 

, n > 1. 



i 

TT 





? 

TT 

? 



TT J 



The data {A„g,Q;„g} satisfy conditions 1-2 of Theorem 1. Let us show that they do 
not satisfy condition 3, and consequently, they cannot be spectral data of L . The 
relations ^{\nq)oinq = , n > , q = l,m for this example can be rewritten in the form 
7(A) = [71(A), 72(A)] , 7i(Aoi) = 7i(Ao2) = 7i(^') = , 72(71^) = , n>l. Clearly, if 
we put 7i(A) = , 72(A) = , we arrive at a contradiction with condition 3. 

Below we investigate condition 3 in some special cases. 

Example 2 (full multiplicities). Let A„i = A„2 = ... = Xnm —'■ A„ for all 

n > . Then rank anq = m , and each of the linear systems ■j{\nq)o:nq = has the 
unique solution 7(A„) = . Wc get the situation similar to the scalar case, because in 
view of asymptotics (3), {A„}„>o can be treated as eigenvalues of some scalar problem. 
Therefore, condition 3 holds automatically. 



We will say that the relations 'y{Xnq)oinq = , q = l,m arc separated for some fixed 
n , if they yield 'yq{\nq) = for all q = l,m . For example, they are separated in the 
case of full multiplicities, or when the matrices anq have a proper diagonal form. 

Example 3. Let the relations 'y{Xnq)ctnq = be separated for all n> tlq . Then each 
component 7g(A) has zeros {\nq\n>no ■ K 7(A) is the function from condition 3, each 
7q(A) cannot have more than no additional zeros (counting with their multiplicities). 
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Otherwise we consider its zeros as the eigenvalues of a scalar problem and prove that 
7,(A) = . 

If 7(A) is entire, and 7(A) = 0(exp(|lm \/A|7r)), |A| — )> 00 , its order is not greater 
than 1/2 . Therefore, by Hadamard's factorization theorem 7q(A) can be presented in 
the form 

7,(A) = (C,o + C,,X + C.^X' + ... + C,,„,A"»)P,(A), P,(A) = J] (l - A 

n>no ^ 

We substitute this factorization into ^{\nq)o!.nq = 0, n<no, q — l,m, and obtain the 
system of linear equations with respect to C^o , Cqi , . . . , Cg^o > g = 1, m . 

More precisely, let Ai , . . . , \n be the first = (no + eigenvalues, and let 
«! , . . . , be the corresponding residue-matrices. For each j = 1, , we choose a 
non-zero column Vj of aj . In case of a group of multiple values among Aj , j = 1,N , 
they have a common matrix aj , and its rank equals their multiplicity, and we choose 
linearly independent columns. Consider N x N matrix P with the columns 

[vjiPi{Xj),VjiXjPi{Xj), VjiX]°Pi{Xj), ...,v jm-fm(Aj), VjraXjPmiXj) 1 . . . , VjmXj -Pm(Aj)], 

j = 1,N . Clearly, that the condition 7(A„g)Q;„g = 0, n <no , q = l,m is equivalent to 
the linear system with the matrix P . Each solution of this system corresponds to 7(A) , 
satisfying condition 3 of Theorem 1. Thus, the condition 3 is fulfilled iff the determinant 
of P is not zero. 



3. Solution of Inverse Problem 1. 



3.1. Let the spectral data A of the boundary value problem L e A{uj) , co E D , he 
given. 
Denote 

D(,, A, = <'^''"f:;'"-'>> = / ^-(t, PM., A) dt. (12) 



We choose an arbitrary model boundary value problem L = L{Q{x), h, H) e A{uj) 
(for example, one can take Q{x) = -u , h = 0„i , H = 0„i). We agree that if a certain 
symbol 7 denotes an object related to L , then the corresponding symbol 7 with tilde 
denotes the analogous object related to L . Put 

m p ma+i—1 p ms+i— 1 p 

Cn ^ ] I Pnq Pnq \ + ^ ^ ^ ^ I Pnq Pnrris I ^ ^ ^ ^ I Pnq Pnms I ^" ^ ^ 1 1 '^^i ^ ^11" 
5=1 s=l q=ms s=l q=ms s=l 

According to Lemmas 1 and 2, 

/ oo \ 1/2 oo 



vn=0 / n=0 



Denote 



XnqO Xnq, Xnql A„g, PnqO Pnijj Pngl Pnij; <^ngO '^ng) '^ngl '^ng' 
V^nqK'^) = Kqi), 'f'nqi{x) = (f{x, Xnqi) , 
Fklj,nqi{x) — (^'klj^i.^^ ^nqij ^klj)j Fklj,nqi{x) — (^'klj^i^J ^nqij ^klj)j 

n,k>0, q, I — 1,171, — 0,1. 
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By the standard way (see [3, Lemma 1.6.2]), using Schwarz's lemma, we get 



Lemma 6. The following estimates are valid for x e [0, tt] , n,k > , r, s — l,m , 
rrir < q < rrir+i , rris <l < rris+i , i, j = 0, 1 : 

\\iPnqi{x)\\ < C, WVnmA^) ix)\\<CCn, 
\\Vnqi{x) - ipnmAx)\\ < C^n, \\Fklj,nqi{x)^ < J^^^^, 



nis+i-l 



— |n-fc|+l' 



l=ms 

\\Fklj,nqi{x) ~ Fklj,nmri{x)\\j \\-Fklj,nmro{x) ~ Fklj,nmrl{x)\\ — ^n-k\+l^ 
{FklO,nqi{x) 

l=ms 

J2 {Fkl0,nmr0{x) ~ FklO,nmrl{x) ~ Fkll,nmro{x) + -^ita,nmrl (^)) 
l=ms 



— |n-fc|+l' 
— |n-fe|+l" 



The analogous estimates are also valid for <fnqi{x) , Fkij,nqi{,x) . 

The lemma similar to the following one has been proved in [7] by the contour integral 
method. 

Lemma 7. The following relations hold 

oo m 

ifi{x, A) ^ (p{x, X) + Y^ ^{(pkio{x)akioD{x, A, Awo) - (pkii{x)akiiD{x, A, Awi)) (13) 

k=0 1=1 



D{x, A, fi) - D{x, A, /i) = ^ ^{D{x, Xkio, n)D{x, A, A^o) - D{x, Xkn, n)D{x, A, Xkn). 

k=0 1=1 

Both series converge absolutely and uniformly with respect to x & [0, tt] and X , /i on 
compact sets. 

Analogously one can obtain the following relation 



oo TO 1 



$(x. A) = $(x. A) + 2^ 2^ }_^{-iy<pkij{x)akij- 



k=0 1=1 j=0 

It follows from Lemma 7 that 



A- A 



klj 



(14) 



<fnqiix) = (fnqiix) + ^ '^i(PkloFklQ,nqiix) - ^kllFkll,nqi{x)) , (15) 



fe=0 (=1 



Frfpu,nqi (•^) Fj^p^^^^qi [ X 



) — ^ 'y^X^m^M^i^)^klQ,nqi{x) — F^^kll{x)Fkll,nqi{x)) (16) 



A;=0 1=1 



for n, r;>0, q,p = l,m, i,uj = 0,1. 
Denote 

^o{x) = Yl (-'^y^Pkiji^WkijVkiji^)^ s(^) = -24(^)- 

{k,i,j)ev 



(17) 



14 



Using (5) and Lemma 6 one can easily check that the series in (17) converges absolutely 
and uniformly on [0,7r] , and the function eQ{x) is absolutely continuous, and the com- 
ponents of £{x) belong to L2(0, tt). 

Lemma 8. The following relations hold 

Q(x) ^Q(x)+e(x), h^h-eo(0), H^H + Soin), (18) 



Proof. Differentiating (13) twice with respect to x and using (12) and (17) we get 

oo m 1 

(f'{x, A) - eo{x)(f{x, A) = ip'{x, A) + ^ ^ "^{-ly <f'^ij{x)a'^ijD{x, A, Xkij), 

k=0 1=1 j=0 
oo m 1 

^"{x, A) = ^"{x, X) + J211 J2(-^yW^HiAx, A, A,,,) 

fc=0 1=1 j=0 

+'^'P'kiji^Wkij'filij{^)'fi{^^ A) + ipkij{x)a'kij{(filij{x)(fi{x, A))']. 

We replace here the second derivatives, using equation (1), and then replace (f{x,X) , 
using (13). This yields 

oo m 1 

Q{x)ip{x, A) = Q{x)(f{x, A) + ^ ^ '^{-iy[(pkijix)a'kij{(flij{x), (p{x, A)) 

fc=0 1=1 3=0 

+'^'P'kij{xWkij'filij{x)'fi{^, A) + cpkij{x)a'kij{^lij{x)(fi{x, A))']. 

CanceUing terms with (p'{x, A) we arrive at Q{x) = Q{x) + e{x) . 
Further, 

oo m 1 

^'{0, x)-{h+ somm = u{ip) + Y.^^^y^i^'^^^ykuDio, x, a«,) = 

fe=0 1=1 j=0 

Since ip{0, X) = Im , <^'(0, X) = h , we obtain h = h — eo{0) . 
Similarly, using (14) one can get 



$'(7r. A) + (i/-£o(7r))$(7r. A) = V{^) + 

k=0 1=1 j=0 



'^''^ X — Xkij 



For j = we have V{(pkio)a'^io = 0^ . For j = 1 

i^Ux), ^x, A))|.=. = V*{^l,,Mn, A) - <^^a(7r)V^(l>). 

Recall that V($) = 0^ , V{^) = 0^ and Q;^^y*((^^^) = 0^ . Consequently, we arrive at 
^'(Tr, A) + (i/-£o(7r))$(7r, A) = 0^ . Together with V{^) = 0^ this yields H = H+ein) . 
□ 

For each fixed x e [0, tt] , the relation (15) can be considered as a system of linear 
equations with respect to (finqi{x) , n > , q — 1,171 , i — 0,1 . But the series in (15) 
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converges only "with brackets". Therefore, it is not convenient to use (15) as a main 
equation of the inverse problem. Below we will transfer (15) to a linear equation in a 
corresponding Banach space of sequences. 



3.2. Denote Xn '■= C7 ^'^^ and Xn = for = . Let V he a set of 

indices u = {n,q,i) , n > , q = l,m , i = 0,1 . For each fixed x G [0, tt] , we define 
the row- vector ip^x) = [ipu{x)]uev and the matrix R{x) = [Rv,u{^)]v,uev , v = {k,l,j) , 
u — {n, q, i) , by the formulae 

ll^nms0{x) = Xn{<fnms0{x) - <fnmA{x)), 1pnmsl{x) = ^nmsl{x), 

'4^nqi{x} Xn{'^nqi{x} ¥'nmsi('^))) 
ms+i-1 

,nmrO 

(x) 

Z=ms 

nis+i-l 

RkmsO,nmrl{x^ Ck FklO,nmrl{x) , 

l=ms 

Rk'rnsO,nqi{x^ Xn^k {,FklQ,nqi^X^ FklQ,nmriixy) , 

l=ms 

Rklj,nmr0{x) = { — iyXn^k{Fklj^nmro{x) — -Ffcij>mrl (^)) ) 

Rklj,nmrl{x) = {-ly CkFklj, mrl{x), \ (19) 

Rklj,nqi{x) ( ly Xn^k{^F}^ij ,fiqi{x^ Ff;ij^njji^i(^xy , 
ms+i-l 

nnirO^X^ FklO,nmrl 

(x) 

l=ms 

—Fkn,nmro{x) + Fkn,nmrl{x)), 

m^+i-l 

R'kmsl,nqiix) — Xn E {,FklO,nqiix) -^feiO,nmT-i('^) -^A;U,nqi('^) ~l~ -^feU,nmT-«('^)) ) 

l=ms 

nis+i-l 

E {FklO,nmrl{x) — Fkii 

,nmrl 

l=ms 

n,k > 0, r,s — l,p, rus <l < m^+i, < q < rrir+i. 
Analogously we define iIj{x) , R{x) by replacing in the previous definitions (Pnqi{x) by 

(finqiix) and Fklj,nqi{x) hj Fklj,nqi{x) . 

We will also use a shorter notation. Consider the row vectors with matrix components 

ipn{x) = [^Pnw{.x), (^„ii(x), (Pn2o{x), iPn2l{x), ■ ■ ■ , <Pnmo{x), Vnm\{x)\, 
l/jnix) = [llJnlo{x),'ll^nll{x),1pn2o{x), 1pn2l{x), . . . , ll^nmo{x) , l/^nmlix)], U > 0, 



and defined analogously 2m x 2m matrices Fi^^{x) , Rk,n{x) , n, /c > , F^; 



klj,nqi 



.X = 



(—iyFf;ij^nqi{x) ■ Then definitions (19) of ij^nqi(x) and Rkij,nqi(x) can be rewritten in 
the form 

^|Jn^^nXn, Rk,n ^ F^^^X^, n,k>0. (20) 

where X„ is a 2m x 2m matrix with components determined from (19). Analogously 
we define <fn{x) , ipn{x) and F^^[x) , Rk,n{x) ■ Now we can rewrite (15) and (16) in 



the form 



<fn^Vn + Y^ ri > 0, 

k=0 



(21) 
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F~ — F~ = \^ F^ (22) 

r),n ri,n / j ^ ri,k k,n \ I 

By virtue of Lemma 6 

\\Rklj,nqi{x)\\, \\Rklj,nqi{x)\\ < |^ _ ^ (^3) 

where C does not depend on x, n, q, i, k, I, j 

Let a„ , u E V , he m x m matrices. Consider the Banach space B of bounded 

sequences a = [a„]„gy with the norm \\a\\B = sup ||a„|| . It follows from (23) that for 

uev 

each fixed a; G [0, vr] , the operators I + R{x) and I — R (here I is the identity operator) , 
acting from B to B , are linear bounded operators. 

Theorem 2. For each fixed x e [0, tt] , the vector ip{x) e B satisfies the equation 

i^ix) ^ iP{x){I + R{x)) (24) 

in Banach space B . Moreover, the operator I + R{x) has a hounded inverse operator, 
i. e. equation (24) is uniquely solvable. 

Proof. Using (20) we get 

(fin = V'n-'^n ^ Fk,n ^ ^kRk,nXn^, 

Substituting these relations into (21), we derive 

OO CO 

i^nX-' = V'n^-' + Yl ^kX^'XuRk,nX-^ = i^nX-^ + iPkKnX;^' ^ > 0- 
fe=0 fe=0 

Multiplying the result by , we arrive at (24). 
Similarly we get from (22) that 

OO 

Rri,n Rri,n — ^ ^ Rr),kRk,n- 
k=0 

This yields R(x) - R{x) - R{x)R{x) = , i.e. {I - R{x)){I + R{x)) = I . Symmetrically, 
one gets (/ + R{x)){I — R{x)) — I . Hence the operator (7 + R{x))~^ exists, and it is a 
linear bounded operator. □ 

Equation (24) is called the main equation of the inverse problem. Solving (24) we find 
the vector ip{x) , and consequently, the functions (finqi{x) by formulae 

Vngi{x) = (PnmA^) + Cn^nqijx) , (25) 

n >0, s — l,p, rus < q < m^+i, i = 0, 1. 

Then we construct the potential Q{x) and the coefficients of the boundary conditions 
h and H via (18). Thus, we obtain the following algorithm for the solution of Inverse 
Problem 1. 
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Algorithm 1. Given the data A. 

(1) Choose L G A{uj) , and calculate ip^x) and R{x). 

(2) Find ip{x) by solving equation (24), and calculate (Pngi{x). 

(3) Construct Q{x) , h and H by (18). 

4. Sufficiency. 

4.1. Let data {^nqiCinq}n>o,q=T^ ^ Satisfying the conditions of Theorem 1 be 
given. Choose L e A{uj) , construct Tp{x) , R{x) , and consider the equation (24). 

Lemma 9. For each fixed x e [0, tt] , the operator I + R{x) , acting from B to 
B , has a bounded inverse operator, and the main equation (24) has a unique solution 

tpix) e B . 

Proof. It Is sufficient to prove that the homogeneous equation 

P{x){I + R{x)) = 0, (26) 

where f3{x) = [f3u{x)]uev , f^uix) are m x m matrices, has only the zero solution. Let 
(3{x) e be a solution of (26), I.e. 

Pnqi{x) + Pklj{x)Rklj,nqi{x) = 0^. 

{k,i,j)ev 

Denote 

lnmsl(x) = /3„m,l('^): 7nm,,o(.«) = lnmsl{x) + ^n/3nmsO{x), 

n > 0, s = l,p, rus < q < m^+i, i = 0,1. 
Then juqiix) satisfy the relations 

oo m 

^nqi{x) + 2^ ^{'yklo{x)FklO,nqi{x) - 'fkllix) Fkn,nqi{x)) = Qm, Tl > 0, (27) 
fc=0 1=1 

and the following estimates arc valid 



\hnqiix)\\ < C{x), n > 0, q = l,m, 

)\\<C{x)U (28) 

s = l,p, rUs < q < rus+i. 



Construct the matrix-functions 7(0;, A), r(x. A) and B[x,\) by the formulas 
lix, Aj = - > }.\lkiQ[x)akio 7 7 7wi(^j«wi r r , (29) 

. — : L A — Akin A — Aun J 



k=0 1=1 



r(x. A) = - > yhkio{x)akio 7 7 lfkii{x)akn r r , (30) 

t^ol^i A-Afc,o A-Afea J 

B(x,X) ^-f*(x,X)r(x,X). 

In view of (12), the matrix-function 7(0:, A) Is entire In A for each fixed x . The 
functions r{x, A) and B{x, A) are meromorphlc In A with simple poles Xnqi ■ According 
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to (29), 7(a;, Xnqi) = 7ngi(a^) . We calculate residues of B{x, A) (for simplicity we assume 

that {Xnqo} r] {Xnql} =0): 

Res B{x, A) = 7*(x, X„qo)^{x, Xnqo)oinqo, Res B{x, A) = 0^. 

A — AngO A — Anql 

Consider the integral 

/iv(2;) = 2^ / ^^^'^^ 
where Fjv = {A: |A| — {N + 1/2)^} . Let us show that for each fixed x e [0, tt] 

lim I]y{x) = Om- 

Indeed, it follows from (12) and (29) that 

CO p ms+i— 1 

-7(a;, A) = ^ ^ ^ ^jkio{x)a'^iQD{x, A, Afc^o) - jkii{x)a'f,i^D{x, A, Afcu)] 

fe=0 s=l (=m3 

oo p rug+i — l 

= 5^X1 [(7fcm.o(2:) - Jkmsl{x)) "W0^(^' ^kio) + lkmsl{xWk\^i^ ^ ^fem.o) 

k=0 s=l l=ms 

rns+i-l 1 

-D{x, A, Afc^,i))+7fcr„,i(a;)(a[.*^-4'^)L'(a;, A, Afc„,i)+7fcr„,i(a;) ^ ^^^^•(^'(a;. A, Aw^) 

l=ms j=0 

mg+i-l 1 

-£>(a;. A, Xkmsj)) + Y^^''i^^^'> ~ lkmsjix)WkijDix, A, Afc^^) 

l=m3 j=0 



By virtue of Lemma 6, (5) and (28), we get 

6 



Ux,X)\\<C{x)eM\r\x)Y \p_k\ + ^ 



, Rep>0. 



Similarly, using (30) we obtain for sufficiently large p* > : 

||r(x. A) II < ^ exp(-|T|x) Y \pJk\ + r ^^P^^^ \p\ ^p*^p^ 

Then 

"TV- 



This estimate yields lim In(x) — Om ■ 

N^oo 

On the other hand, calculating the integral iNi^) by the residue theorem, we arrive 



at 

oo m 
fe=0 q=l 
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Since am = ali^ > , we get 

lkio(^hkio{x)akio = 0^, 



7(a;, Xkio)oikio = ^m, k > 0, 1 = 1, m. 
Since 7(2;, A) is entire in A , and 

7(3;, A) = 0(exp(|r|a;)) 

for each fixed x G [0, tt] , according to condition 3 of Theorem 1, we get 7(a;, A) = 0„i . 
Therefore '~inqi{,x) = 0^ for all n > 0, g = l,m, i = 0,1, i.e. the homogeneous 
equation (26) has only the zero solution. □ 

4 -2. Further, we provide the general strategy of the proof of sufficiency in Theorem 1. 
The proofs of Lemmas 10-12 are similar to ones described in [3, Sec. 1.6.2]. 
Let — ['^u{x)\uev be the solution of the main equation (24). 



Lemma 10. For n>0, g = l,m, i = 0, 1, the following relations hold 

V'„,.(a;) eCi[0,7r], < C(n + l)^ z/ = 0, 1 a;e[0,7r], 

where 

Construct matrix-functions (Pnqi{x) by formulae (25). By virtue of Lemma 10, we 
have 

yi%ix)\\<Cin + ir, u = 0,l, 

ynqiix) - >pnqiix)\\ < 0^7]^, y'nqiix) - >Pngiix)\\ < CT] , g = 1,771, (31) 
\\Vnms0{x) - (PnmA^)\\^ W^nqiix) - (Pnmsi{x)\\ < C^n, S = l,p, TUg < q < m^+i. 

Further, we construct the matrix-functions </7(x. A) and $(x. A) by the formulas 



ip{x, A) = (^(x. A) - ^ {-iyipkij{x)a[ 



ik,i,j)ev ^ - ^''^ 



$(a:,A) = $(a:,A)- J] [-1]' ^kiA^)a'u, 'Y,^ ^^ 

and the boundary vahic problem L{Q{x), h, H) via (18). Clearly, ip{x,Xnqi) = (Pnqi{x) ■ 
Using estimates (31) one can show that the components of EqIx) are absolutely con- 
tinuous and the components of e{x) belong to L2(0,7r) . Consequently, we get 



Lemma 11. Qjk{x) G L2(0,7r) , j,k = l,m . 
Lemma 12. The following relations hold 

i^nqiix) = Xnqi^nqiix), £(p{x, A) = X(p{x, A), £^{x, A) = A$(x, A), 
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(^(0, A) = Im, ^'{0, A) - h, U{^) = Im, V{^) = 0, 



In order to finish the proof of Theorem 1 it remains to show that the given data 
{Xnq,Oinq} coincide with the spectral data of the constructed boundary value problem 
L{Q, h, H) . In view of Lemma 12, the matrix-function A) is the Weyl solution of 
L . Let us get the representation for the Weyl matrix: 

{k,l,j)eV 



oo m 



A;=0 



A — Xkii A — Xkii 



Using the equality (see [4]) 

oo m 



Ol'i. 



A — Afc/i 

fe=0 i=\ 



«fe«0 



we arrive at 

oo m 
fc=0 «=1 '^"^ 

Consequently. {A^/o} are simple poles of the Weyl matrix M(A) , and {aao} are residues 
at tlie poles. Note that the multiplicities of the eigenvalues coincide with the numbers 
of equal values among {Afe^o} > because they both coincide with the ranks of {afejo} ■ 
Theorem 1 is proved. □ 
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